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1. Introduction 

Over decades, the model of a localized magnetic impurity in a non-magnetic metal has 
been one of the major challenges in many-particle theory. Anderson pQ proposed a model 
of a localized impurity interacting with a host of free electrons through hybridization: 

H A = ^2 e k c k,T C k,r + e d n d + ^2 ( M fc4 c fc,r + M£cl T d T J + Un dA n dd 

k,r k,r 

On the impurity site, a Coulomb repulsion U is allowed. The scattering at the impurity 
is assumed to be isotropic and therefore one- dimensional. In the limit \Mk\ 2 /U <C 1, 
a localized moment forms, which is demonstrated by a canonical transformation of the 
Anderson model, resulting in the Kondo model [2] with an impurity operator 

Hi = 2 J S ^2c\ T (T TiT iC k i y(Ti , (1) 

k,k' 

where cr = (a x ,a y ,a z ) denotes the Pauli matrices. The Kondo model describes a 
free host, interacting weakly with a localized magnetic moment via antiferromagnetic 
XXX spin exchange with an amplitude J. "Weak interaction" means that at high 
temperatures, the coupling is negligible and the impurity spin shows Curie-Weiss 
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behavior. The lesson to be learned from this limit is that a localized moment occurs if 
singly occupied sites are energetically favourable and hybridization only leads to virtual 
double or zero occupation. The model (£[} served as starting point for Kondo J2J, who 
performed a perturbational calculation of the scattering amplitude between host and 
impurity up to third order in J. He discovered a In Tk/T contribution to the electrical 
resistivity. Tk is the crossover temperature which indicates the limit of perturbation 
theory: A divergence occurs for T ~ Tk- 

The method which overcomes the failure of perturbation theory is scaling. By 
the implementation of his numerical renormalization group, Wilson achieved a non- 
perturbative calculation of the impurity contribution to the magnetic susceptibility \ 
and the specific heat C at T <C T^. By assuming a linear dispersion in the conduction 
band he discovered Fermi-liquid-like behavior of the impurity for T <C Tk- In the other 
extreme, at T ^> Tk, Wilson confirmed the asymptotic expansion in In Tk/T discovered 
by perturbation theory techniques. 

Andrei [Sj and Wiegmann [Zj obtained the spectrum exactly by the Bethe Ansatz 
(BA). The linearized energy- momentum relation turned out to be crucial for the 
application of the BA. Thermodynamic equilibrium response functions were calculated 
in the following by employing thermodynamic BA (TBA) techniques, jHj. The impurity 
contribution to the free energy is encoded in a set of infinitely many coupled NLIE. These 
contain the whole information about thermodynamic equilibrium functions. Especially, 
the asymptotic high-temperature expansion due to Kondo and the Wilson ratios are 
encoded therein. The low temperature Fermi-liquid-like behavior was confirmed in the 
cited works. However, the high-temperature asymptotic expansion was not performed 
as far as in Wilson's approach. 

We develop a lattice path integral representation of the free energy of a one- 
dimensional Anderson-like impurity model in a correlated host. This model can be 
viewed as a lattice-regularized version of the Anderson model in the continuum. As a 
special case, the Kondo model ist obtained in a certain scaling limit. The host is based 
on a four-dimensional representation of the Lie superalgebra gl(2|l). The corresponding 
four states per lattice site are zero, single (with spin up or down) and double occupation. 
The impurity degrees of freedom are described by a three-dimensional representation of 
gl(2|l). Double occupation on the impurity site is excluded from the beginning. The 
parameters of the model can be tuned such that on one hand, particle exchange between 
the impurity site and adjacent host sites is eliminated by a canonical transformation 
and on the other hand zero occupation is energetically suppressed. The same parameter 
tuning makes vanish correlations in the host. These conditions fulfilled, a localized 
moment in a free host occurs (Kondo limit). 

In order to regularize the continuous Kondo model, it is quite natural to choose 
the superalgebra gl(2|l). Its even subalgebra is u(l)CS>su(2), encoding charge and spin 
degrees of freedom, respectively. Spin-charge separation occurs in one dimension for 
interacting electron systems and the impurity is supposed to possess exclusively spin 
degrees of freedom. Indeed, the scaling limit reduces gl(2|l) to one of its subalgebras, 
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su(2), in the impurity space. Then the excitation spectrum contains only spin degrees 
of freedom on the impurity site. 

The model proposed in this work allows for the Kondo limit as one special case. We 
calculate the free energy exactly in the general case, that is an Anderson-like impurity in 
an interacting host and perform the Kondo limit afterwards. Thus our results are farther 
reaching than the known non-perturbative treatments of the Kondo model j3J|H]. The 
free energy of the host and of the impurity are given by eigenvalues of distinct quantum 
transfer matrices and can therefore be separated. In the Kondo limit, Wilson's results 
are confirmed with high accuracy. The general Anderson-like case will be investigated 
elsewhere [TT)] . 

This article is organized as follows. In the next section, we derive the Hamiltonian 
by QISM. The third section deals with the calculation of the free energy. In each of 
these sections, the Kondo limit is treated explicitly. Section four contains the derivation 
of Wilson's results in the framework of our path integral approach. A conclusion and 
an outlook form the last section. 

In all what follows we set fee = 1, and g\i# = 1, where fee is Boltzmann's constant, 
g is the gyromagnetic factor and \xb is the Bohr magneton. An index % (h) denotes 
quantities pertaining to the impurity (host). 



2. The impurity model 

Let V^ d > be the module giving rise to the ci-dimensional irrep of gl(2|l), d — 3, 4. A 
grading is assigned to the basis vectors through the parity function p, 

d = 4 : p[l] = p[4] = ; p[2] = p[3] = 1 
d = 3 : p[l]=p[2] = 0; p[3] = 1 . 



(2) 



The matrices R 



(d,d') 



'•J 



GEndfV^ ® V/ d M satisfy the graded Yang-Baxter-Equation 



(YBE), 



Rtf\u) 



Pa 

P' a' 
(d",d) 



R\ 



(d",d'), 



oca 



V — U] 



a', j3' 



R\ 



(d",d') 



J a" ,P" 

&'a' 



a ,7 



a" ,7' 



tiffin) 



(p[a']+p[a"])p{f3'] 



13" a" 



(3) 



Summation over doubly occurring indices is implied in the foregoing equation and in all 
what follows. 

Explicit expressions of the R matrices are given in the following, 

' \p[a}p[b}„b ^ a) 



R^ 3 \u) = (u + (-l)fWP[6] e 6 g, e «) 

R^\u) = u I (u + | + 1 + (-ipMpMe* Et) 



R^\ 



u 



2a 



u — a 



-Pi 



2a + 2 
u + a + 1 



(4) 
(5) 

(6) 
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e a (-^a) are the nine three- (four-) dimensional generators of gl(2|l), obeying 

[ e b' e d]± := e t> e d ~ (~ l)( p ^ +p ^( p ^ +p ^ e c d e^ 

= Vd e t ~ {-l)^ [a]+p[h ^ [c]+p[d ^5te a d , (7) 

and the same for the E£. Pi, P3 are projectors from <S> onto gl(2|l) modules 
with highest weights (0, 0\2a) and (—1, — l|2a+2) respectively. They are given explicitly 
in For a matrix representation of e£, choose the basis 

|T) = (1,0,0), |2) = (0,1,0), |3) = (0,0,1). 

Then e£ := \b)(a\ is the usual matrix representation of projectors in three dimensional 
space. 

As to the E£, we choose a basis in V^ A \ 

|1) = (1, 0, 0, 0) , |2) = (0, 1, 0, 0) , |3> = (0, 0, 1, 0) , |4> = (0, 0, 0, 1) . 

We call the projectors associated with these states m£ := |&)(a|, a,b — 1,2,3,4. One 
verifies that the set [TT] 

P 1 1 = -|3)(3|-|4)(4|, P 2 2 = -|2)(2|-|4)(4|, 

El = a\l)(l\ + (a + 1)(|2)(2| + |3)(3|) + (a + 2)|4)(4|, 

E^\2)(3\, P 2 1 = |3)(2|, (8) 

El = v^|l)(2| + v / ^TT|3)(4|, El = y/a\2){l\ + v / ^TT|4)(3|, 

E\ = -Va|l><3| + v / ^TT|2)(4|, El = -\/a|3)(l| + V / ^TT|4)(2| 

satisfies eq. (J7J). In the sequel, the real parameter a is restricted to a > 0. 



Consider the set of matrices R defined by 



R (u, 



a,0 
-,.6 



4)P[<5](P[7]+PM) R( d '' d )(- U 



7,/3 



(9) 



The permutation of the indices means that creators and annihilators are exchanged in 
the auxiliary space of R. These P-matrices satisfy 



R [u) 


Pa 


R {d "> d '\v) 










a'ri" 



R (d",d) { 



v — u\ 



a,f3 
a', j3' 



R K 



R(d",d) (v _ 



u , 



a', 13' 



a", (3" 
— (d,d') 



«) 



"I) 



(p[a]+p[a'])p\j3'] 

_l\(j>W]+p[a"])p\J3'} _ ^ 



The R- matrices can be translated into a graphical language. Straight lines denote 
the 4-dimensional space, wavy lines symbolize three-dimensional space, 

P 



[R (4A) (u-v)] 



y.8 = a 
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The two lines symbolize the two spaces intertwined by Rij. Each line carries a 
direction indicated by an arrow; both the vertical and horizontal lines carry spectral 
parameters. The argument of R is given by the difference between the right and the 
left "incoming" parameters. The replacement R — > R means flipping the arrow on the 
vertical bond. 

The YBE eq. (J3J) in graphical language reads: 





"Other" YBEs are obtained by flipping arrows (that means replacing R — > R) and/or 
substituting straight by wavy lines (that is, changing the dimension in one of the spaces). 
"Unitarity" is a further property of the i?-matrices. 



R«< d 'h 



u 



/3,7 
5,a 



-u 



a,S 



-y',/3' 



8 P P 



11^ 



and the same for R. The unitarity property fixes normalization constants of the R- 
matrices. In the following, we will speak of "normalized" /^-matrices when they satisfy 
eq. (fTTj) ; non- normalized /^-matrices differ from those by constant pre- factors, but still 
fulfill the YBE. A direct verification of eq. (jllj) for d = d' — 3 (d = d' = 4) is done by 
using the projection properties 

e a b e c d = 5 a d e c b , Pj P k = 5 jtk P k , (12) 

or by using the YBE with standard initial conditions. Furthermore, for d = 3, d 1 = 4 in 
eq. (fTTj). one should employ 

4)p[0](p[a]+p[<5]) — _( a _|_ 2)(— l) p Hp[*l£; Q . 



In order to construct a lattice model, one embeds m b e b „ into End 



such that e b a acts non-trivially only on the lattice site 0. Therefore consider the graded 
tensor product of two operators v, w. 

V a b ® s Wd = (-l)(pl»]+Plb])p[d] v £ w c ? 

where v,w stand for e, m. The operator of unity in three- (four-) dimensional space is 
J 3 = e c c (J 4 = m£). Following [T2], define 



(_ 1 )(p[o]+p[6])EiLiP[cfc] e « g, m ci g 

® B (L-j) 



(8> m 



<-L 



with j 



— ^_i)(p[ a ]+pM)Efe^ + iP[cfc]j 3 
. , L. Then 

M = (-l)W^W)W^M)[m fc ]S[eo]JJ . 



(13a) 
(136) 
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Analogous relations hold between rrij, m^. 

Principally, at this point one could derive the Hamiltonian. However, it is more 
convenient to find a fermionic representation of the R matrices in order to use the more 
familiar language of fermionic field operators c| ■, c T j, acting on the spin directions 
t =1, 1 and on the lattice site j. This is done by employing the technique of Gohmann 



[T2l ITS] , which consists in identifying the 
fermionic operators. 

The entries [Xj]^ of the matrix 



in 



Jib'' 



[ e o\b with certain combinations of 



X j 



( n n n^ 



j 

-it 

-Cil(l-^'T) 



(1 - n n )n^ 



-n i; )(l-n 



„t t 



(14) 



it) / 



satisfy projection and commutation properties formally identical to eqs. (|13a|) 
with grading p[l] = p[A] = 0, p[2] = p[3] = 1 in accordance with eq. © for d = 4. This 
is the only constraint on so that we identify [Xj]^ = [rrij]^. The whole set 

reads in fermionic language: 

[ E £ = a + 2- (n n + n iT ) 
\ E j\\ = n H ~ 1 

ml = -V® n n c ji - 

[Ejt = V® n m\ ~ v / « + 1 (1 - n n )c n [Ejf 2 = y/E n n c^ - s/a + 1 (1 - n n )c ] n 



1 

1(1 



1 — n 



ml 
ml 
mi 
mi 



n 3l 



an tf c ji ~ V" + 1 (i — n jd> 



The even sub-algebras are manifest: £3 is the u(l)-generator, and E\\ are the su(2) 
generators. The fermionization of the three-dimensional e£ is done with the matrix Y, 
resulting from X, eq. (fT4^) by deleting the first row and column, 



Y 



' n d)l (l-n dA ) 
d[d\ 

V 



n d,i) n d,\ 



-d\{l-n dA ) 
(1 -n d>l )d\ 



\ 



(15) 



(! -nd,i)(l -n dA ) J 



-di(l-n dA ) (l-n da )d T 

We slightly modified the notation, replacing c\ c, n by d\ d, n d . Horizontal and vertical 
bars separate fermionic and bosonic sectors. The boxes on the diagonal of Y contain the 
generators of su(2), u(l). Set [e } b a = Y^, such that eqs. (|13aj) . (j!36|) hold with grading 
{1,1,0}. 

The monodromy matrices 

T(u) = Rtt\u)RttM ■ ■ ■ R a A f{u)R a f{u + m ) (16) 

T W=R a ,L {- u )Ra,L-l{- u ) ■■■R a ,l {- U ) R a,0 (~ U + m 0) 

consist of sequences of R matrices, multiplied in (horizontal) auxiliary space. Note the 
shift by iuo on the zeroth lattice site, where the dimension of the (vertical) quantum 
space is reduced by one. This site shall be denoted as "impurity site". The shift is done 
by iu G C, for reasons which will become clear later. Graphically, T(u) is depicted as 
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i 

5' 



L 1 

The super-trace is called transfer matrix 

t(u) = str a T(u) , r(u) = str a T(w) (17) 
In [tt] (u) = In [rr] (0) + tt [r- 1 (0)r'(w) + r- 1 (0)r , (u)] u=o +0(u 2 ) . (18) 



=: const. if 



In the last line, the Hamiltonian was defined as the logarithmic derivative of the two 
transfer matrices at zero spectral parameter. By scaling tt, one is free to multiply the 
Hamiltonian by a constant factor. 

Before evaluating eq. (jlHjl . let us shortly comment on the case of a homogeneous 
model without impurity. We denote the corresponding quantities with a subscript 
h. This model has been extensively studied in ^TJ Ej • Assuming periodic boundary 
conditions 1 = L + 1, r h (0) (j h (0)) is the right (left) shift operator, and 

lnr^O) =iP = -In 77,(0) , (19) 

where P is the generator of translations to the right. The derivative with respect to u in 
eq. (|T%|) yields a sum of L terms, each one corresponding to Rj j +1 (0), j = 1, . . . , L. For 
the ease of notation, let us follow the graphical depiction of jT^]. The following figure 
shows the jth term of [ r /T 1 (0) r /l( M )] u= o : 



h 



jj+i 




M0) 



j+i j 



The vertex with a dot denotes i?'- +1 (0). For tt = 0, the vertices decouple and taking 
the trace over a row yields the right shift operator 7^(0) . Thus 

L 
3=1 

-ff/j is scaled by D(a + 1), D is a bandwidth parameter whose significance will become 

>(4,4) 



clear later. Using explicit expressions for -R^j+i from ^Tj one confirms the expression 
for i 

%j+i 



for /ijj+i given in [T7j 



a 



a 
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D E( c Wi,r + c] +liT c ilT )e-3^.^+^ 

T 

+ U ( n jJ n j,l + n j+i,l n j+i,l) + t P (c]+i,T c ]+i,l c i.T c ja + c )A,l c i+ l ^ c i+ l >i) 
+ D(n j + n j+ i) -2D , (21) 
a + 1 



U= — = t p , e - " 



r = — r . 



O: 



One can show by a canonical transformation ^7] that the limit a — » leads to the 
t — J-model. On the other hand, for a>l, the leading order are free fermions: 

L L 



1 + 2a Jli ' f + n 3 +1 ^ 



H h /D =2j2nj-2-J2I2 ( c lr c i+i,r + ct +ljT c jV 

2 L 1 L / \ 

+~ Z^ n J.t n i,i + - E ( c ]+i,T c ]+i,l c i.T c ia + c i,T c L^+ 1 .T £ ^+ 1 .V + ( D /« 3/2 ) • ( 22 ) 

Since in this work, we aim at realizing a free host, our interest is in the limit a>l. 

Due to the insertion of R^'o , H receives an impurity contribution H^. It can be 
derived graphically. First observe that eq. (JT§j) no longer holds; but due to unitarity 
JTTJ), one still has r(0) = r _1 (0). Then r(0) is depicted as: 

i 

£' 

c 



V \ 



1 L 

In comparison with the free host, the changes induced by the impurity stemming from 
lnr'(O) correspond to the graphs: 

L 



i 

s 



1 



1 L 



A vertex with a dot symbolizes the derivative with respect to the spectral parameter. 
The first term, 



p(3,4)/ • \ e>1*>> 



,(4,3) 



p-1 (4,3)/. x 
H 10 [Mo) 



d(4,3) 
-^1,0 



UmoJ 



couples the impurity to the left neighboring site. The second term 

d(3,4) 
^1,0 



-1U 



h L ,l d(4,3)/- , \ _ z-,-l(4,3)/ , h L> i (4,3) 



wtoj 



D(a + l)" l ' u " i ' u v "" u/ D(a + i; 

is a three site coupling. Analogous terms, with L and 1 interchanged, are provided by 
lnr'(O). The inverse matrix [i?( 4,3 )] is found by eq. (jlljl. In analogy to the host 
Hamiltonian, the spectral parameter u is scaled by D(a + 1), where -D is a bandwidth 
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parameter. Using the fermionization technique, we calculate Hi. It contains terms 
which are symmetric and others which are antisymmetric under (L <-» 1) exchange. In 
the combined thermodynamic and continuum limit, the latter do not interact with the 
impurity [TH], so we do not consider them in the following. Then 

Hi = -2DJ a {\ - n d ) {a + F hL ) V d 



- DJ a ^/aV d ^ [ d H c L,r + Ci, T ) - d T {S Lr + c\ . 

T 

+ DJ a ^ n d,r \ Cl, t Ci,t + 4,r C L,r) V d 

T 

+ DJ a 4dr (c^c 1)T + c\ >f c L , T ) P d + O (D Ja/a 1 / 2 ) 



(23) 



where the projector Vd '■= 1 — ^d,t^d,i projects onto non-doubly occupied states on the 
impurity site. Furthermore, 

F l)L = 2-2h- c} iT c L)T - c f u c La - 4 )T ci )T - 4 a ci a 
2a 

J « = 2 , 2 >0^o :=«o/2 (24) 
with the most interesting range of the coupling constant 



a 



2 < J < 



a 



(25) 

In jAppendix A[ an alternative fermionization is used, resulting in essentially the same 
Hamiltonian. 

Finally, one includes external fields /i, h, by 



Hp 



h 



^ rij + n d V d 



(26) 



Eqs. f|22|) . (J23J), (I26j) define the entire Hamiltonian of the impurity model in the limit of 
an asymptotically free host. In |Appendix B it is shown that Hh + Hi displays gl(2|l) 
symmetry. H ex breaks this symmetry, but preserves integrability. 

The Kondo limit can be performed by a canonical transformation which eliminates 
transitions between single and zero occupation of the impurity site. It is conveniently 
done in Fourier space: 



1 

7? E C k,r' 



s/L 



l k t 



2 x ^ 



;q,q' 



Q ( ,Q\ t t 

2^ cos -cos {q+-j c q+QjT c q>T c^,_ Qf c q%f 



2(cos£; + 1) -fi/D 



(27) 
(28) 
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Hi = {2DJ a (n d -l)(a + F 1)L )-fin d }V d 

+ J a {aD) l ' 2 V d Y,{M k d\c k , T + Mlc\ T d T )V d 



+ J a ^ N k,k' ^,r4 T c fc / iT + d\drc\ -c k < )T V d (29) 



T,k,k' 

M k = - ^ (1 + e ifc ) , N Kk , = j (e" ifc ' + e ifc 

Here I = L/D is the constant length of the chain, and D~ l plays the role of a lattice 
constant. The canonical transformation is generated by an operator A, which yields a 
transformed Hamiltonian H e ff = exp(A)H exp(— A) not containing any hybridization 
between impurity and host in the leading order 0(J Q ). One verifies that 

H^E (M k d\c Kr - M*4 r d T )P d , (30) 

, Q — Cfc 

where := 2J a an d — njD has been defined. H e ff contains terms O(aJ^). Given the 
restriction (|25|) . these terms can be neglected. After the transformation, the excitation 
spectrum of the impurity site contains only the contribution for single occupation. The 
contribution for non-occupation in the impurity operator is energetically suppressed 
in the strong coupling limit: In the language of renormalization theory, Vq drives the 
impurity Hamiltonian to a strong coupling fixed point at low temperatures. This will 
be demonstrated in the next section. Thus we do not consider the contribution from 
zero occupation to the Hamiltonian in the ongoing. 

To perform the scaling limit, the fermionic spectrum is linearized around 
incommensurate Fermi points ±hp avoiding Umklapp scattering. The linearization gives 
rise to right (left) moving particles R (L). To avoid divergences due to the unbounded 
linear spectrum, operator products are normal ordered: 

where v e {R, L}. The continuous description is achieved by introducing field operators 

1 L 

J - _ V P ifc - vt 

* n=l 
L 



C n,u,r^^ D 

n=l 



1 



e 



VlJo 

where x = n/D, vpl T (x) = limo^ vDc^ VT , q k = k ■ D. ijj are now fermionic 
field operators with|^, )T (x), ip\, T , (x') j = 5(x — x')5 u y5 TtT i. Again normal ordering is 
imposed, 

e— >(J 
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where (• • -)o is the expectation value in the ground state. Let us summarize the external 
fields again in an operator H ex . Then 



Hh / n r . sin kp j. . . d . , 

— = / o v y(cosk F + 1) :n v : - b v y\ v + w— —\ty\ v {x)—ty T y{x)\ 

** ,.,./_ L 



Hi = 2cosk F J / 5(x) y~" j \o~ v yn d)T V d :n v ^{x): + d\d f :ipl )f (x) ip v >y{x): 



dx 



dx , (31a) 
(316) 



H ex = / -fj, [n{x) + 5(x)n d V d ] + - [5(x)(n d j - n d>i ) + (n T (x) - n;(z))] dx . (31c) 
Jo * 

The occupation number operators are n TjU = :tpl u ip T y., n v = J2 T n T ^, n T = ^2 u n T>l/ . 

As far as the terms ()31 ojl . (|316J) are concerned, one may pass to a Weyl basis by 

the canonical transformation 

Interaction terms in the host are non-local in the 0±(x); however, as will be shown 
in the next section these are accounted for by a redefinition of the Fermi velocity Vp, 
smkF = vp — > vf = vf(1 + 0(l/a)). The Weyl basis demonstrates that the impurity 
couples only with one of the two host channels. The Hamiltonian density thus reads: 



iv F --(j)]. T {x)—(j) v Jx)\ + D(cosk F + l):nJx): - Dt u 
dx 



(32) 



Hi = AJcoskp ^2 $( x ) '■ ( ^+A x ) ( t ) +A x )'- n d,T'P d + :(jJ + ^ T (x)(j) +i f(x):d\d T 



H ex = -(j, [n{x] + 5(x)n d ] + - [6(x)(n d ,i - n dA ) + n T (x) - n^x)] . 

The fermionic operators of the impurity can be expressed in terms of spin operators 
with index i, 

<? z i = n dA - n dA , af = d\d t , a~ = d]d T . 

Then one directly recognizes that the impurity operator is su(2)-symmetric and can be 
completed to the XXX-exchange operator, 

Hi = 2J5(x) 22 ■(t)+A x ) cr r,T'(p+,T'(x):cri + 2J5(x)n d V d :n + (x): , (33) 

t,t' 

where cr = (a x , a y , a z ) T and 2 cos kpJ a ='■ J is defined. Eq. fl3^|) constitutes the 
isotropic Kondo model. 



3. Calculation of the free energy 
Taking account of eq. (fTSj). 



e pHh = lim [T h (u N )r h (u N )] N/2 , u N 

N—*oo 



N/2 
N/2 -/3H e 



pD(a+l) 
N 



e~ pH = lim [f{u N ) r(u N )] N/2 e 

N-^oo 
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e ~(3H ex _ J~J e -/3[h/2(n JiT -n ja )-/inj] e -f3[h/2(n dA -n dd )-fiY, T n d ,r] _. e ~PT, 



f =1 h e 



-PH e 



J'=l 



The even integer N is referred to as Trotter number and is the height of the fictitious 
underlying square lattice, see fig. (JIJ. The impurity contribution to the free energy is 



fi = — lim lim — 

L->oo N^oo f3 



lntr 



[t{u n ) t(u n )} 



N/2 -fiH ex 



lntr 



[f h (u N )r h (u N )} N/2 e-^ 



where H ex ^ = "52j=i hex,j- The crucial idea in calculating fi is to exchange tr and str in 
the expression 

N/2 

tr { [r(u N ) t(u n )] n/2 e-P H °* } = tr e~ pH " fj str 

k=l 



-UN 



x R afo{- u N + mo)i?i 4 2 ^ lL (u7v) ■ • • R i ^ ) _ ll {u N )R i ^ ) _ i0 {u N + iuo) 



This leads to 

L 

str J 



x 



AT/2 



tr^-e 



(4,4) 



"AT J 



fc=i 

7V/2 



tr e-^n<;o(-^ + i%)^ (4 ' 3) 



Q2fe- 



i0 (ujv + m ) 



(Q), 



k=l 

N/2 



--: str 



T<%„) 



tr, e 



tr n e~ 



-I3h e 



if 



fc=l 

AT/2 

-/3J? e x,i TT"p (3,4) i 

0,a 2fc ' 



(Q), 



if 

fc=i 



^ + b)i?£ 4) -(«jv + iu) =: tr T, 



(Q), 



(34) 



(35) 



Eqs. (|34|) . (J35|) define the Quantum Transfer Matrix (QTM) of the host and 7$ of 
the impurity, respectively. Note that the host matrix is independent of the lattice site 
j. Each QTM is the trace over the auxiliary space of a Quantum Monodromy Matrix 
r p(.Q) m The auxiliary space of is four-dimensional, of t\ three-dimensional. Fig. Q 
depicts this "rotation" from auxiliary space into quantum space. Due to eqs. @, (fTUj) . 



r^(v),4?\v')\ =0, (36) 

where the symbolical indices u, v' may take values h, i. The auxiliary spectral parameter 
is essential for the diagonalization of the u^- are inhomogeneities with alternating 
signs. Especially, eq. ()36|) holds for v ^ v'\ The impurity and host QTM's share the 
same set of eigenvectors. The largest eigenvalue of tP is separated by a gap from the 
rest of the spectrum for any N. The eigenstate |$ max ) leading to the largest eigenvalue 
A™ ax (w) of Ti(v) also leads to the largest eigenvalue A™ ax (w) of (u). Although 
interesting, this is not essential: The dominant eigenstate l^max) of the host matrix 
determines the "correct" eigenvalue of the impurity matrix 
l 



In 



str ( 7^(0)) r™(«o) 



M, 
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Figure 1. Classical lattice representing the free energy of the impurity model. L is 
the physical, N the Trotter direction. The dimension in the impurity quantum space 
(wavy line) is reduced by one. The spectral parameter on the straight vertical lines 
is and on the wavy vertical line — iw . For the further analysis it is convenient to 
introduce the auxiliary spectral parameter — iv on the vertical lines. Crosses stand for 
twisted boundary conditions, induced by external fields h and /i. 



Ill 



111 



;-l) p[max] (A™ ax (0)) A,- 



E 

fc^max 



(-l)'M (a«(0)) A?>(«o) 



(A^(0)) L Ar x K) 



+ 



fc^max 



.-j\p[fc] 



A?Vo) 



Af(0) 

A™ ax (0) J Af iax (w ) 



(37a) 



(376) 



which is an asymptotical expansion for large L. Generally, the eigenstate of the kth 
largest eigenvalue of does not lead to the kth largest eigenvalue of t^. So with 
respect to t\ , k does not label the eigenvalues according to their order. The supertrace 
requires to include the parity of the projector on the eigenstate k. Note that p[max] = 0. 
In JH] it is argued that the two limits N — > oo, L — > oo are interchangeable. Then 
the thermodynamic limit L — > oo in eq. (j376|) is carried out just by keeping the largest 
eigenvalues A™| x . 

One concludes that the impurity and host contribution to the free energy per lattice 
site are given by 



fi 

fh 



lim — In A' 
lim — In A 



max 
h 



(0) 



(38) 
(39) 



Eqs. (|3*H|). (jSni) summarize the enormous advantage of considering the QTM: The 
calculation of the free energy is reduced to the evaluation of a single eigenvalue. 
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fh has already been calculated in this approach [20*1 I21|. the result is given below. 
The diagonalization of is done by applying techniques of the nested Algebraic Bethe 
Ansatz (NABA) jT^j, yielding for the eigenvalue Aj of the non-normalized QTM t\ . 

Ai(v) = X.(v) + X + (v) + X (v) (40) 
X_(v) = q -^ V ^ <P + (v + ia/2)0_(v - ia/2 - i) e ^ +h ^ 

X + (v) = g+ ^~ " ia/2U + (v + ia/2 + i)eP^~ h M 

= g-fo + i)q + (v - i) Mv _ ia/2)(f) _ {v + ia/2) 
q_(v)q+(v) 



q+{v) 



M M 

]J(v - v s ) , q-(v) =H(v- v k ) , Mv) = (v ± m) N ' 2 . 
j=i k=i 

An external magnetic field h and a chemical potential fi have been introduced. The roots 
- or particle solutions - {v^} are determined by the analyticity of the eigenvalue: 

K(vj) = q-( v ) (j>-(v-ia/2) <j> + (v + i + ia/2) ^ {fJt _ h/2) 
\{ v j) Q-(v + i) 4> + (v -ia/2) cj)_(v + ia/2) 

^-{vk) = q+(v) (j>+(v + ia/2) <j>-(v -i-ia/2) ^ +h/2) 
X (vk) q+(v-i) 4>+(v - ia/2) + ia/2) 



-1 (41a) 
-1 . (416) 



v=v k 



These Bethe- Ansatz equations are M + M many nonlinear coupled algebraic equations 
for the unknown roots. Using analyticity properties, we represent the eigenvalue by 
a finite set of non- linear integral equations (NLIE). Within this approach, the Trotter 
limit N — > oo is carried out analytically. Consider the following combinations of \± t o, 
eq. (jUJ): 

1 X+(v) / , X (v) 



b(v) ■ A_(«) V A+(«0 



q+(v - i) <t>-{v - ia/2) <t>+{v + i + ia/2) fih 



q-(v + i) <p + {v + ia/2) <p_{v - i - ia/2) 
x (i + + <l>+(v- ia/2) (j>-(v + ia/2) c _^„ fe / 2) 

q+(v)\ q-(v) (f>-(v- ia/2) (f) + (v + i + ia/2) 



V 

q_ {v) 



<j)-(v - ia/2) <p + {v + i + ia/2) 



q^>(v)e-^ , (42) 
4> + {v + ia/2)0_(i; — i — ia/2) q_{v + i) 

where qL 1 := Y[f=i M+M f u ~ and vf^ are called hole solutions. The polynomial 

q_ has been identified by reasons of analyticity: The zeroes of numerator and 
denominator cancel as far as the particle solutions are concerned, the hole solutions 
rest as zeroes of the numerator. The polynomials in the denominator are the same 
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as the 0-terms in Ao/A + . Along the same reasoning (or simply by taking the complex 
conjugate and h — > —h), we find another function b: 

1 A-(u) A , A (u) 



b(u) ' A+(v) V A_(v) 



— ia/2)0 + (t> + i + ia/2) g + (t> — i) 
where := f|^ M+M fi> — uj^V A third function c is introduced 



:^Y>K\ (43) 



1 Mv) Hv) 



c(v)' A + (*;)A» 

<p + {v - ia/2) 0_(u + ia/2) e" 2 ^ 



t^H. (44) 



0_(w - ia/2) 0+(w + i + ia/2) <p + (v + ia/2) 4>_{v - ia/2 - i) 
Consider 03 := 1 + b, «8 := 1 + b, £ := 1 + c, 

»(„) = j^bWUv) = , q -^\ h) Mv)e-^ 
A -{ v ) q + [v — 1) q_ '{v) 

*(«) = Y^rbCt,) A^) = - A^) e -^+ fe / 2 > 

C ( B ) = 1 c(t> ) = ggM^ U 2 ^ 

V7 b(v)b(i;) - ia/2) + ia/2) A^v) 

In Appendix B[ the unknown functions q±, q± \ A$ are expressed through the auxiliary 
functions by means of analyticity arguments for the largest eigenvalue. The result is a 
closed set of NLIE: 

\nb(v) = <f) b {v + i8) - [k b *\nV$}{y + 2i5) - [k b * In €](v + i5) + (3(n + h/2) (45a) 
In b(v) = (j^{v - i8) - [% * In *B](v - 2i5) - [% * In C](v -iS)+/3(fi- h/2) (45b) 
In c{v) = (f) t (v) - [k b * ln»](v + iS) - [% * ln»](u - id) - [k c * In €](v) + 2/5/i (45c) 
The contributions of the impurity and host to the free energy are given by: 

-Pfi= - In c(u ) - 2pD(a + 1) J a + 2/3 fi (46) 
-13 f h = 77(0) + [C * In 05] (0) + [C * In ®] (0) + [(C + C) * In <£] (0) (47) 
The inhomogeneities are 

Mv) = - - f ^-^v +1)2 /o > h = ft (48) 

(f + ia/2)(t> — ia/2 — 1) 
0c = 0t, + <k ■ (49) 

The convolutions [/ * g](x) := J"^° f(x — y)g(y)dy involve local kernels: 

1 

2irv (v — i) 



k b {v) — — - — , k b — k b , k c — k b + k b 



D(3(a+1) ' /w y v 2 + (a + l) 2 ' 
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Following the treatment of the Hamiltonian in the preceding section, we want to 
perform an asymptotic expansion of the free energy in the limit a ^> 1. The 
essential observation from the study of the Hamiltonian is that after the canonical 
transformation, excitations of the impurity stem exclusively from transitions between 
singly occupied states. Furthermore, it was argued that uq has to be scaled such that 
zero occupation is energetically suppressed. First consider the case where Uq is held 
fixed. By scaling v — > av, the algebraically decaying kernels shrink to 5-functions. 
The leading contribution thus solely stems from the driving terms <j) b ^ c . Within this 
approximation, the auxiliary functions can be calculated explicitly: 

„ , \ a(av) -. . a(av) . . a(av)a(av) ,_ n . 

b(av) = -, — r , b(av) = -. — r , c(av) = - -. — r -. — - (50) 

l + a(av) l + a(av) v ' 1 + a(av) + a(av) 

a(v) = exp [4> b (v + i/2) + /3(jji + h/2)} , a(v) = exp [^(v - i/2) + p(p - h/2)\ . 

The quantity 6 in eqs. ()45aj) - ()45c|) has been chosen 5 = 1/2 such that a, a are real 
valued functions. Define lima^oo^a + 1) J a =: Jo. From the expression of fi, eq. (jjB)) it 
then follows that the free energy is that of an uncoupled impurity, 

lim fi(T,h) = -Tin \{aa)-\u ) + a _1 (u ) + a" 1 ^)] +2DJ -2fi 

a— >oo 

= — Tin [eP 2DJo + e^ +h ' 2) + e^~ h ^} . (51) 

The free energy reflects the expected result of a free impurity with three possible states, 
namely empty and singly occupied with up or down spin. A more detailed analysis 
of the case where Uq is held fixed will be given in a forthcoming publication JU]. We 
now demonstrate that the two latter states dominate if uq is scaled appropriately with 
a, D, fi, such that a crossover temperature emerges, below which a strong coupling fixed 
point is reached. 

For low temperatures, the auxiliary function c exhibits a sharp crossover from c 1 
to c ^> 1 in regions around "Fermi points" ±A C defined by 



-0 c (A e )«2/3/i, K^±a^-- A . (52) 

Set h \i. The influence of h on the Fermi points is neglected, since it enters 
quadratically. The more common parameterization is 

v 1 k ,_ nN 
— = - tan — , (bo) 
a 2 2 ' K ' 

where k is the wave-vector used in the Fourier representation of the Hamiltonian in the 

preceding section (not the Fourier variable conjugate to v). At v — A c , eq. is 

equivalent to 

2T>(cos£; F + l) = /i , (54) 

which defines fcj? at constant // at D 3> T, such that fi = 6f (/•* in turn is related to the 
particle number N by kp = ttN/(2L). Then eq. (|54|) yields a relation between \x and 
N, independent of T - this demonstrates that the formally grand-canonical description 
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is effectively canonical, because D ^> T). For the analysis of the NLIE, consider 

k e & * In €(v) — k c b * In c(v ) + / k c t,(v — w) In — — r du> . (55) 

y^i<A C ' w 

To evaluate the second integrand, note that asymptotically, 

£ = 4(l + a)(l + a) = o(-^ (56) 



c aa \bb 

11 c(u) [0, |u| > A c " 1 j 

Consider the case |u| < A c in eq. (|57|) . The first term is exponentially small, the second 
term dominates. It is inserted into the second integrand on the rhs of eq. ()55j) which 
shall be considered as a next-leading correction compared to the driving terms in eqs. 
()45aj) - (j45cj) . In view of eq. (|46|). the most interesting range of the spectral parameter is 
v ~ v . Set 

|«o| = K'o/2| > A c ~ oi , (58) 
as indicated in (|25j). Then one proceeds with eq. (|53j) 



k c b * ln « fc cb *lnc(u) + fc cb (v) / ln ^-^dw; + 0(l/v 4 ) . 

As an estimate for the second term on the rhs, one uses the leading term of eq. (|57j) : 
ln^f = [-0 C (^) - + O (exp[-/3£>]) 



In -^dw « 4/foL> 

w|<A, W 



2A C uA c ^. , . 
2 arctan — - - + 0(1 a) 
a olD 



-:(3k>0. (59) 



The sub-leading order 0(1/ a) is neglected in the following, which is justified rigorously 
below. The above defined quantity k is a monotonously decreasing function of n/D. 
Choose A c > such that k > 0. Summarizing, 

kb jC * InC = kb >c * \nc + /3K,kb, c , (60) 

so that ln € can be eliminated in eqs. (j45aj) - (j45cj) . The Fourier transforms of (j>b, (fa, 4>c, 
eqs. (ggj, (gHD are: 

f -a/2fe i; > 



^(fc) = -PD(a + l] 



e -(a/2+l)fc j fc >Q 

e Q / 2fc , fc < 



-a/2k (i _i_ p -fc\ £• > f) 

=^(fc)+^(fc) = -^D(a + l)< l /2 w),t^ j ' ,"n (61) 



e a/2fc (i + e fc ) , k < 



Inserting eq. ([HOjl gives, using eqs. (j45a|) - (l45c|) : 
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b(k) = - - 1 J k + T-^-nu (® - 8~) (62b) 
v ; 1 + e" fc 1 + el fc l 7 v ' 

t(k) = - —^9* + — ^(8" - 8) (62c) 
1 + e* 1 + 

Note that in the limit A c — > 0, i.e. yU — > AD, the resulting equations are trivially solved; 
In b = 0h = - In b, and 05 = 1 + exp(/3/i), 8 = 1 + exp(-/3/i), 

fi = const. + In (e^ 2 + e"^ 2 ) , (63) 

which is the free energy of a free, uncoupled spin. This situation corresponds to a 
completely filled band. 

The NLIE eq. ()62ajl - ()62c|) are transformed back to direct space, 

\nb(v) = -2nf3K$(v + i5) + (3h/2 + [jfe * In 8] (u) - [fc * In 8] (v + 2iS) (64a) 

In b(v) = 2it(3k$(v - iS) - (3h/2 + [k * ln8](w) — [k* ln<B](v - 2i5) (646) 

In c(u) = -/3-D(a + 1) 9 - = - - + 

f z + ar/4 

+ [$*In«8](v-i£) - [$*ln8](t; + io") . (64c) 
The driving term and integration kernel read: 



2 sinh 71 v 



1 e -|fc|/2 



k(v) = — / — — e lKv dk . 

K ' 2ix 2 cosh k/2 

Choose 5 = 1/2 and scale v by 1/ir. Since $ decays exponentially, it is possible to 

absorb k and vq in a new additive constant. Substitute 

v = x — ln(27Tft) (65) 

and remember that k scales with aD and therefore may be arbitrarily large. All 
parameters can be combined in the free energy to a new constant Tk, 

- \n(2nK) - irvo/2 =: - \nT K , T K = 27r«;e™ o/2 . (66) 

The range of \v \ has been identified in eq. (|58j). we take v = — \v \. The shift (|63j) 
scales the driving term <$>(v/ir — i/2): 

- 2P<irK$(v/ir - i/2) = - (67) 

coshf 



^x— ln(27TK) _|_ x+\n(2-KK) 

= - f3e x . 

In the second line, eq. ()65j) has been employed. At this point it is clear that sub-leading 
orders in eq. (jSUjl can safely be neglected. Furthermore, a factor ft can be absorbed by 
shifting x — > x — ln/3: 

lnb(x)= -e x + (3h/2+ [fc*ln8](x) - [A; * In 8](x + ivr - ie) (68a) 

lnb(x) = - e x - /5/i/ 2 + [fc * ln8](x) — [k * ln8](x - ivr + ie) (686) 

1 Z" 00 [ln88](a;) , , , 

-Pfi = 7T / —? -^^dx . (68c) 



2tt 



cosh ^a; + In J 
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As far as the host is concerned, eqations (|5U|) are inserted into equation (J47)) : 

- hm (3f h = v (0) + ±- [ Ac/a ^—\ni(l + a(v))(l + a(v))}dv. 
a^oo 2vr J Ac/a v 2 + 1/4 

Observe the relation between the elementary excitation energy e(v) and the momentum 
k(v) as functions of the spectral parameter v (cf. eq. (|53JI ). 

e(v) = ^-k(v) . (69) 

Then one obtains the energy-momentum relation e(k): 

e(v) = -= r- — > k(v) = 2arctan2t> 

k 

e(k) =4 cos 2 - = 2 cos A; + 2 (70) 

The function k(v) in the first line is given by eq. (J69|) . From eqs. (J59j) . (}66j) . Z} gets 
arbitrarily large, leading to a linear dispersion in the host. Then the free energy density 
fh of the host is given by: 

rp pkpD 

fh = _ Hm f- / In [1 + e-^M-W)] [i + e -/^M+V2-M)] dg ( 71 ) 

/3>1 27T J^ kpD 

Jl := fi - 2D I ( 1 + - ) (coskp + 1) - -(cosA; F + l) 2 
\\ a J a 

Vf '■— 2 ( 1 H (cos fci? + 1) ] sin fc^ , 

\ a a J 

where q = k ■ D. Note that interactions in the host of order 0(l/a) can be absorbed 

into a redefinition of vp = 2sinkp, resulting in effectively free fermions. The leading 

orders of the specific heat and magnetic susceptibility are 

vr 2 2 7T 2 

C h {T)=T--^=:T-p h (72) 

3 TTVp 3 

Xh(T) = i =: , (73) 

4 7rUp 4 

where is the density of states in the host.f 

4. Calculation of high- and low-temperature scales 

We demonstrate the advantage of our novel approach by a direct calculation of high- 
and low-temperature scales and comparison with Wilson's results jl]. 

Wilson found for h,T <C the ratio of the specific heat Ci to the magnetic 
susceptibility Xi is universal: 

1 7rT 



2nT K 7 y ' 3T K 
_ X i(T) C h (T) 

Rw -cm~x~m- 2 - (74) 

f The host in our model, eq. I|32(l . contains two channels. Therefore ph is enhanced by a factor 2. The 
impurity couples only to one of the two channels. 
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We confirm these results in our approach by using dilogarithmic identities j22] to extract 
the lowest order of the free energy for low fields and temperatures, 

T 2 7T 1 h 2 

lim f(T,h) = --=- . (75) 

From eq. (JZEJ, 

lim C(T) = — £ lim = 7^ • (76) 

This Fermi liquid behavior is to be compared with the host, eqs. (|72|). ()73|) . Define the 
coefficient of the linear T-dependence of Ch by 5^. Then R^, is calculated as 

:= lim ^ = 2 . (77) 

For high T and ft, = 0, the impurity susceptibility asymptotically reaches the Curie- Weiss 
limit: 

1 



X(T) 



1 lnx , _ 3 x 
O (x d ) 



x=ln T/Tjf 



X 2x 2 

Tk is defined such that contributions O ^ln~ 2 T/Tk) do not occur, 

T K = 2n£T K , £ = 0.1032 ± 0.0005 , (78) 

where the numerical value of £ was calculated by Wilson The calculation of the 
numbers R w , £ are benchmarks for any non-perturbative solution of the Kondo model 
covering the entire temperature axis. In [Hj, a perturbative expansion of the free energy 
by Andrei gave an analytical expression of £, £ = 0.102676 . . .. This number has not 
been obtained in the framework of the exact TBA solution yet. However, by analyzing 
the NLIE eqs. ()68a|) - fl68cj) we are now able to give an accurate numerical value of £, 
which is summarized in the rest of this section. 

In the high-temperature regime T ^> Tk, the asymptotic behavior of the auxiliary 
functions for x ~ — IhT/Tk <C gives the main contribution according to eq. ()68c|) . 
The convolutions in eq. ()68aj) are evaluated asymptotically at h — 0: 

7r 2 

Re ln«B(x < 0) = In 2 + 



Ax 3 

1 1 In \x 

2 4 x 8x 2 ' X' 
_ 1 / , 1 In |s| , 40+1 



Re d ph \n*{x «0) = - + — - ^ + ^ (79) 



Re^ln«( a; «0)- i ^--—-r— - 

The coefficient of the x _2 -decay in eq. (|T9|) is determined numerically; we find 
(f) = 0.04707±2-10~ 7 . From this one gets £ = exp [-40 - 1/4] /2vr = 0.102678±2- 10" 6 , 
agreeing nicely both with Wilson's and Andrei's results. Details of the calculations will 
be published in a forthcoming paper [TO] . 
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5. Conclusion and outlook 



An Anderson-like impurity Hamiltonian on a one-dimensional lattice has been obtained 
as the logarithmic derivative of a gl(2|l) symmetric transfer matrix. The free energy 
for the lattice model was calculated exactly from a closed set of finitely many NLIE. 
As a special case, this impurity model allows for the Kondo limit of a localized 
magnetic impurity in an interacting host of electrons including the free fermion case. 
Mathematically, the Kondo limit constitutes the reduction of gl(2|l) symmetry to su(2) 
symmetry on the impurity site. 

We expect that the lattice regularization presented in this work can be generalized 
to the Kondo model with anisotropic exchange by using the quantum super-algebra 
U 9 gl(2|l). Furthermore, it is possible [2H] to find representations of gl(2|l) such that 
one of its subalgebras is the (2S + l)-dimensional irrep of su(2). This allows for a lattice 
path integral approach to multichannel spin-S* Anderson-like models. In a forthcoming 
paper [10], we will give the corresponding closed set of NLIE in the Kondo limit. These 
can be obtained by an argument of symmetry, and allow for the calculation of Wilson 
ratios in the general multichannel spin-5 case. After having found a lattice path integral 
approach in the simplest case (spin-1/2, single channel), the next question addresses 
two-point correlation functions like (crfa^), r = 1,...,L. Since the spectrum and 
the corresponding eigenstates of Tf®' are known, those quantities can be calculated 
in principle, by generalizing methods developed in to graded models. In view of 
progress in describing the screening cloud around the impurity at T T% by methods of 
conformal field theory |2lH |2*4~] complementary understanding from the exact solution is 
highly desirable. Finally, the Anderson impurity model on the lattice, including charge 
fluctuations, can be treated exactly by solving eqs. ()45aj) - (j45c|) . a question currently 
under investigation. 

We acknowledge very useful discussions with F. Gohmann. 
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Appendix A.l. Alternative fermionization 

The matrices X and Y defined in (1 14 J) . (fT5j) provide one possible fermionic representation 
of the [mj] b a , [eo] b a . This representation is not unique: In this appendix, we use slightly 
modified matrices X, Y. However, the fermionic representation of the Hamiltonian is 
essentially the same. Consider 

/ (! - n il)( l -"it) (i-^'lta'T c ii(l 



X j 



t 



~ c jl c tf 



n 3l0- ~ n jh 



n 



(AA) 
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V 


-d\n dA 







Again one may identify [mj]? 



Y 



since eqs. ()i3fo 



are fulfilled. 



The matrices X, X are related by a particle-hole transformation, so that fermionization 
of Hh yields the same expression as in eq. (|2~T|) . The fermionic representation of the 
impurity operator is denoted by Hf. 

Hi = -fm d + ^ (n d)T - n dj i) - 2DJ a n dA n di (a + F 1>L ^ 

- DJ a y/a ^2 Uti 1 ~ n d,-r){cL,T + ci )T ) - d T (l - 7tyr)(4, )T + c\ T ) 

T 

- DJ a ^ n d,r \ 2n lt r + 2n L ^ - c ] L -c l:T - c\ -c Lt ^ 

T 

+ DJ a J2 dU- (2c{ >r c L>T + 2c{ T c 1:T - c{ >r c lyT - c{ 7 c LjT ) + O (DJ a /a 1/2 ) , 

T 

where 

2F ljL = 2(m + n L ) - c f ljT c LiT - c\ A c LA - 4 >T d, T - c[ 4 c u . 

The canonical transformation generated by A, eq. ()30|) . is also applicable to Hi, and 
the doubly occupied state is energetically suppressed by scaling w . In the Kondo limit, 
Hi and Hi are seen to be identical. 

Appendix A. 2. gl(2\l)- symmetry 

From the YBE (jSJ), the direct product of two monodromy matrices is intertwined by an 
i?-matrix, 



',13" 



= [R^\v - u)]% [T^\v) T^\u)] a X„ . (A . 2) 

The invariance of r(u) with respect to gl(2|l) is shown by expanding eq. (jA.2|) in the 
limit v —* oo, keeping only terms 0(1), 0(l/v). 



1 + I (f + 1 + (-1)^3 



v \2 
>( 3 > 4 )/„,m?(3,4) 



o U 



(A.3) 



Ra7(v)R a T 1 (v)...R a f(v + iu 



- 1 + £ {e [f + 1 + t e J« w + [e a ) b a K 

V \v* 



O 



(A.4) 



T^ 4 ' 4 ^(ti) = T(m) is defined in eq. (FTB|) . 
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Eqs. (|A.3|) . (|A.4J) are inserted into eq. (jA.2|) with d" = 3; d — d' = 4, while keeping 
the full T(it). The constant terms on both sides are identical. In order 0(l/v), one gets 

( _ 1)p[/ 3»](PW +P K]) [T(n)] J ;I y« + [T{u) f pii [(-iJPl^]^ 

= (-l) p[/3](p[a]+p[a " ]) ^[^(^)]^ + [(-l) p[aMb] eX]^ • 

Set /? = /?", multiply with (—I) 13 and sum over /?. The second terms on each side are 
identical. The first terms give the commutator of the transfer matrix t(u) with W: 

t(u) := [T(u)g , [r(«), W£] = . (A.5) 

Dropping constants in W, eq. ()A.5|) states that: 

L 

r(u),J2m a b + [eo] a b =0. 
j'=i 

Thus r commutes with all global gl(2|l) symmetry operators. In a very similar way, one 
starts with eq. (fTUj) to show 

n^x>,]: + M =o, 

3=1 

where t{u) is defined in eq. (TTTj) . Consequently, the Hamiltonian, defined by eq. (JTEJ), 
is gl(2|l)-symmetric. 



Appendix B. Derivation of NLIE 

The unknown functions in eqs. ([40 )1 - (|44 [) are q + , g_, ,q^ ,q+\ where the indices 
pertain to two different sets of particle and hole solutions. Incidentally for the largest 
eigenvalue, the index denotes the part of the complex plane where these functions have 
zeroes: If a q— function carries an index + (— ), it has zeroes in the upper (lower) half 
plane. From numerical studies for finite N, we know that in the largest eigenvalue case, 
the "particle solutions" obey Im[f j] > , Im[?) fc ] < Vj, k. Analogously, the "hole 
solutions" are distributed in the complex plane as Im[t^. ] < 1, Im[u£ ] > 1 , Vj, k. The 
particle and hole solutions are distributed discretely and accumulate at certain points, 
prohibiting a formulation in terms of densities. 

In the following, the largest eigenvalue case is studied. Consider the logarithmic 
derivative of these auxiliary functions, so that constant terms vanish. Since we know 
the analyticity properties of all functions in the complex w-plane, we can calculate their 
Fourier-transforms, 

r°° chi 
f = J [In/Wfe- 1 ^. (B.l) 

The integration contour is taken along the real axis. This is allowed as long as 
\ot/2\ > \un\, which certainly is the case for and a sufficiently large. / vanishes 
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for fc < (fc > 0) for f(v) analytic in C + (C~). Thus it is convenient to calculate the 
Fourier transforms separately for k < 0, k > 0. For the moment, concentrate on k < 0. 

_ g(jfc) = _ e^ 2+1 ^> k ^_(k) + e k q + (k) (B.2) 

-t(k) = - e ka/2 4>4k) - e k q + + qf (B.3) 

-i(k) = e fca/2 (0+(fc) - 0_(fc)) - e (Q/2+1)fc 0_(fc) + Ai(k) (B.4) 

«B(fc) = -e^AO+A^fc) (B.5) 

t(fc) = q+{k) - q+\k) + Ai(k) (B.6) 

C(k) = - e ka / 2 Mk) + (fc) - ki{k) . (B.7) 

The essential observation is that in eqs. ()B.2j) - (jB.7j) . there appear the three unknowns, 
namely q + , q+ and A i; and the three auxiliary functions b, b and c. Add eqs. flB.6|) . 
(IB.7|) and combine this sum with eq. ()B.2|) : 

b(k) = e {a/2+1)k (4>4k) - + (fc)) - e k (<5(k) + £(fc)) . (B.8) 

Combine eqs. (jB.5|) with ()B.2|) and these two with eq. (|B.8[) . An expression for A, 
results, 

ki(k) = »(*) + e k (t(k) + €(k)) + e (a/2+1)fc <Mfc) , 
which is inserted into eq. (|B.4|) : 

c(fc) = -e ka/2 (4> + (k) - i_(k)) + e (Q/2+1)fc (0_(A;) - + (fc)) - <8(fc) - e fc (^3(fc) + £(fc)) . 

Finally, eqs. ()B.5|) and ()B.7|) give q+ (k), which is inserted into eq. (|B.3|) . 

t(k) = e fcQ / 2 (0_(/c) - $ + (k)) - (£(fc) + »(*)) . (B.9) 

The case fc > is obtained by exchanging b, b, switching fc — > —k in the exponential 
terms and replacing 0_ *-> <p + . The result is summarized: 

e- ka / 2 {<f) + {k) - }_(k)) - (3*(Jfe) + £(fc)) fc > 
b(fc) = <( -(^(fc) + £(fc)) fc = 

e («/2+i)*(0_(A;) - <^ + (fc)) - e fe (^5(fc) + C(fc)) fc < 

e -K2+i)fc(0 + (A;) - 0_(jfe)) - e - fc (Q3(fc) + <£(fc)) fc > 
b(fc) = <( -(»(*) + C(fc)) fc = 

e W2(0_(jfe) _ - ($8(fc) + fc < 

( e -(a./2+i)* + e -te/2) (0 + (A;) - 0_(fc)) - ^(jfe) - e - fc (<B(fc) + £(fc)) fc > 
c(fc) = { -%(k) - (<8(Jfe) + C(fc)) fc = 

( e (a/2+l)fc + e W2) _ _ _ e fc (^(fc) + £(fc)) fc < 

Application of the inverse Fourier transform and integration leads to a system of non- 
linear integral equations. 

lnb(v) =0jf°(u + W) - [fc**ln»](v + 2W) - [fc 6 * In + iS) + /3(p + h/2) (B.10) 
lnb» = - i5) - [%*ln«8](v - 2i8) - [fc F * ]n€](v - iS) + f3(^ - h/2) (B.ll) 

lnc(u) = 4>[ N) (v) - [fc b * ln»](i; + W) - [fc ¥ *lnQ5](u - i5) - [fc c * In £](«) + 2/3 fi ,(B.12) 
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where the convolution 

lf*g]{x) ■= 

is done with local kernels: 
1 



f{x - y)g(y)dy 



(B.13) 



k b (v) 



k- b ( v ) = k b ( v y 



k c (v) = k b (v) + k%(v) 



27ru(i;-i)' " 0K "' ' '^ y " J " yv "' ' 2tt(v 2 + 1) ' 

In order to achieve convergence, the equation for In b (lnb), eq. (jB.10|) (eq. (|B.11|) ). is 
taken for v + iS, (v — i8). 

The constant terms are integration constants derived from the asymptotic behavior 
of the auxiliary functions for large \v\. 

a aa 



lim b = — , lim b = , lim c 

\v\— >oo 1 -\- a, \v\—*oo 1 -\- a, \v\— >00 

a = e P(»+h/2) - = e /3(»-h/2) _ 



1 + a + a 



(B.14) 
(B.15) 



The inhomogeneities are 



<p[ N \v)= In 



^ + (u + if) 0_ fc-if -i) 
J_ (« + if) + («-if-i) 



The thermodynamic limit iV — >• oo can be carried out leading to eqs. (|45ajl - (l45cjl . 
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